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The cataclysmic variable SDSS J1507+52: An eclipsing 
period bouncer in the Galactic halo 
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^p^' SDSS J1507+52 is an eclipsing cataclysmic variable consisting of a cool, non- 

Q I radially pulsating white dwarf and an unusually small sub-stellar secondary. The sys- 

V-^ , tem has a high space velocity and a very short orbital period of about 67 minutes, 

c/5 ' well below the usual minimum period for CVs. To explain the existence of this pecu- 

I liar system, two theories have been proposed. One suggests that SDSS J1507+52 was 

formed from a detached white-dwarf/brown-dwarf binary. The other theory proposes 
that the system is a member of the Galactic halo-population. 
^ \ Here, we present ultraviolet spectroscopy of SDSS J1507+52 obtained with the 

■ Hubble Space Telescope with the aim of distinguishing between these two theories. The 

00 ' UV flux of the system is dominated by emission from the accreting white dwarf. Fits to 

model stellar atmospheres yield physical parameter estimates of Toff = 14200 ± 500 K, 
log5 = 8.2 ±0.3, usini 180±20 kms"! and [Fe/H] = -1.2±0.2. These fits suggest 
■r^lj- I a distance towards SDSS J1507+52 of d = 250± 50 pc. The quoted uncertainties 

. include systematic errors associated with the adopted fitting windows and interstellar 

reddening. 

Assuming that there is no contribution to the UV flux from a hot, optically thick 

boundary layer, we find a Teg- much higher than previously estimated from eclipses 

analysis. The strongly sub-solar metallicity we infer for SDSS J1507+52 is consistent 

, with that of halo stars at the same space velocity. We therefore conclude that SDSS 

H ' J1507+52 is a member of the Galactic halo. 

Cd ■ 
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1 INTRODUCTION 

Cataclysmic Variables (CVs) are binary systems where a 
late-type main-sequence donor star transfers mass onto a 
primary white dwarf (WD). According to standard evolu- 
tionary theory, loss of angular momentum drives CVs to 
initially evolve from longer to shorter orbital periods. AH 
systems are predicted to reach a minimum orbital period 
roughly when the donor stops its hydrogen burning and be- 
comes degenerat e. The observed minim um period is found 
at « 83 minutes ICansicke et al.l l|2009l ). At about the point 
when the donor reaches this evolutionary state, the thermal 
time scale becomes longer than the mass-transfer timescale, 
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and the donor cannot shrink fast enough in response to con- 
tinued mass loss. As a result, the brown-dwarf donor will 
expand, and the system will move towards longer orbital 
periods. 

Standard evolutionary theory also predicts that about 
70% of all CVs should have pa ssed their minimum pe- 
riod and have sub-stellar donors l|Kolblll993l : iHoweh et al.l 
Il997f ). However, few such systems have been successfully 
confirmed, and until recently, almost no CVs containing 
donors with masses below the hydrogen-burning limit, had 
been found. Systems that have passed their minimum period 
(post-bounce systems) should be very faint and are there- 
fore difficult to find, in general. Nevertheless, in the course 
of the Sloan Digital Sky Survey (SDSS), a few post-bounce 
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CVs were discovered (jLittlefair et al among them the 

ecHpsing system SDSS J1507+52 (hereafter J1507). 

This particular system was quickly recognised to be 
odd due to its sho rt orbital period of about 67 minutes 
ISzkodvet al.ll2005l '). which is well below the minimum pe- 
riod characteristic of normal CVs. A few other systems are 
occasionally found to have periods below the usual minimum 
period, but, in general, these systems show evidence for ab- 
normally hot and bright donor stars in their optical spec- 
tra, suggesting that their secondaries are nuclear-evolved 
objects. This is not the case in J1507, whose secondary is 
not visible at all in o ptical spectroscopy (|Szkodv et al.|[2005l : 
iLittlefair et al.|[2007l V Together with its short orbital period, 
this indicates a system with a faint disc, a low accretion rate 
and a relatively unevolved, low-mass donor. 

ILittlefair et al.l l|2007h performed eclipse analysis of 
J1507 and found a donor mass of 0.056 ± 0.001 Mq, clearly 
below the hydrogen-burning limit. However, because of its 
anomalously short orbital period, the system is not consis- 
tent w ith the standard m ass-period relation for C Vs l|Kniggel 
I2OO6I ). In line with this, ILittlefair et~all (j2007l ) found the 
donor radius to be smaller than predicted by standard 
CV evolution sequences. They therefore suggested that the 
donor might be unusually young, i.e. that J1507 might be 
a CV that formed recently from a previously detached WD 
- brown dwarf close binary system. A young brown dwarf 
has a higher density, and therefore a smaller radius, since 
the donor has not yet had a chance to expand in response to 
mass loss. This would explain the short orbital period found 
in J 1507, since according to the period-density relation, a 
higher density and subsequently a smaller radius, implies a 
shorter orbital period. The list of photometrically inf erred 
system parameters provided by ILittlefair et al.l l|2007h also 
includes the effective temperature of the white dwarf, which 
they estimate to b e Tea = 11000 ± 500 K . 

Concurrently, iPatterson et al.l (|2008l ) showed that the 
system has an unusually high space velocity, similar to the 
velocities of stars in the Galactic halo. A typical star of 
1 M0 in the Galactic disc has a space velocity below 50 
kms~^, while J1507 has a velocity of about 167 kms~^. If 
J1507 is a member of the Galactic halo, its donor will be 
a Population II object with low metallicity. Due to their 
lower atmospheric opacity, such objects have significantly 
smaller radii than their solar metallicity Population I coun- 
terparts. Membership of the Galactic halo would therefore 
also account for J1507's small donor radius and short or- 
bital period. However, Population II stars only represent 
about 0.5 % of the stars in the S olar neighbourhood, w hich 
would make J1507 a rare svstem. [Patterson et al.l (|2008l ) es- 
timated a slightly higher effective temperature of the WD in 
J1507, Tcft = 11500 ± 700 K. They also found the system to 
exhibit multi-periodic variability on time-scales of minutes, 
which they interpreted as non-radial pulsations originating 
from the primary white dwarf. 

Both hypothesis imply that J1507 is an interesting and 
important system that can shed light on poorly understood 
aspects of CV evolution. However, since the models are very 
different, it is clearly important to determine which - if ei- 
ther - of them is correct. As noted above, eclipse-based es- 
timates for the effective temperature of the WD suggested 
Tefi ~ 10000 K - 12000 K. At these temperatures, the ultra- 
violet (UV) spectrum of a WD accreting metal-rich material 
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Figure 1. The Gaussian smoothed and flux-calibrated light curve 
obtained by COS/HST, constructed from the monochromatic flux 
averaged over the region 1425 A - 1900 A. Bottom panel shows 
the phase-folded light curve in the vicinity of the eclipse. 

would be expected to show strong absor ption features due to 
Fe II and Fe HI (c.f. iHorne et al.lll994l i. We have therefore 
obtained far- and near-UV spectroscopy of J1507 with COS 
and STIS onboard the Hubble Space Telescope, with the aim 
of measuring the metallicity of the system and establishing 
its status as either a Pop I (disk) or a Pop II (halo) object. 



2 OBSERVATIONS AND REDUCTION 

Far- and near-UV observations of J1507 were performed in 
February 2010, using the HST. The observations were ob- 
tained in the so-called time-tag mode. In this mode, each 
photon event is tagged with a corresponding event time, 
along with their wavelength and spatial position on the de- 
tector. This allows us to rebin the data over any time interval 
and wavelength region. 

Far-UV observations were carried out with the Cosmic 
Origins Spectrograph (COS), for a total exposure time of 
about 3 hours, covering 3 eclipses in three successive HST 
orbits. The maximum time resolution of the COS TIME- 
TAG mode is 32 ms. The G140L grating was used, which 
gave a spectral range of 1230 A - 2378 A and a spectral res- 
olution of 0.5 A. 

Near-UV observations were obtained by the Space Tele- 
scope Imaging Spectrograph (STIS), for a total exposure 
time of about 1.4 hours, at a maximum time resolution 125 
/is. The G230L grating was used, resulting in a spectral 
range of 1650 A - 3150 A and a spectral resolution of 3.16 A. 
The data set consists of 3 sub-exposures (of length, 1600 s, 
1918 s and 1708 s). No near-UV data was obtained during 
the eclipses. Also, gaps in the data are cased by interruptions 
coming from Earth occultation. 

All data were reduced and calibrated using the Pyraf 
package Stsdas, provided by STScI (the Space Telescope 
Science Institutejj. 

^ Stsdas and Pyraf are products of the Space Telescope Science 
Institute, which is operated by AURA for NASA. 
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Figure 2. The two top panels show the mean out-of-ecUpse spectrum for both the far (left) and near (right) UV regions. A model 
constructed from parameters obtained from the far U V-region is plotted in red (see Table [TJ . A model corresponding to values from 
Littlefair et al. (2007), refi=11000 K and logg=8.5, is plotted in blue. Underneath each wavelength region is a plot showing the residuals. 
The dips at (1460, 1530, 1600 and 1670) A coincide with regions of bad pixels. The l-cr flux errors are plotted in green. 



3 ANALYSIS AND RESULTS 

Figure [T] shows the Gaussian smoothed and flux-cahbrated 
hght curve from the far-UV data. This light curve was con- 
structed by dividing the larger exposures into one-second 
bins, and by average the fluxes over the wavelength range, 
1425 A - 1900 A. The light curve shows sharp, square- 
shaped eclipses, close to zero flux at mid-eclipse, and a very 
flat out-of-eclipse continuum, strongly suggesting that the 
UV flux is dominated by emission from (or very close to) 
the accreting WD. The folded light curve constructed from 
the combined data for all three HST orbits is presented in 
the bottom panel. All further analysis was carried out on 
the unsmoothed data, excluding the eclipses. Ground-based 
optical observations were performed at the same time as 
the HST observations, yielding a magnitude for J1507 of 
g = 18.44±0.02. This is consistent with SDSS, showing that 
the system was in quiescence during the HST observations. 

The two top panels of Figure [2] show the mean out- 
of eclipse spectrum (in black), for both the near- and far- 
UV. The reduction was performed separately for the COS 
and STIS data, and the mean spectrum for each region 
were joined at 1900 A without overlap for a total range of 
1268 A - 3150 A. The COS data was rebinned slightly, to 
0.2 A/pix, to increase the S/N and provide sufficient sam- 
pling of the G140L grating's spectral resolution. Different 
dispersions were retained for the the far- and near-UV re- 
gions (0.2 A/pix and 1.58 A/pix, respectively). The overall 
UV spectrum for J1507 is relatively sparse in spectral lines. 
We identify CII at 1335, Sill at 1526.7 Hell at 1640.5, AlH 
at 1671 and MgH at 2800 A. The strong emission feature 
near 1300 A is almost certainly due to geocoronal emission 
in the O I 1304 A line that could not be cleanly removed by 
the pipeline background subtraction. 

3.1 Spectral Modelling 

A model grid spanning the four key atmospheric param- 
eters, effective temperature (Teff), surface gravity (logy), 
metallicity {[Fe/H]) and rotational velocity (usini), was 
constructed in order to fit the far- and near-UV spectrum. 
Atmospheric structures were calculated with Tlusty, with 



the spectral synthesis bein g done with Synspec l|Hubenvl 
ll988l : lHubenv fc Lanj[l995l ). The grid consisted of; 12500 K 
^ Tefi ^ 20000 K in steps of 500 K, 7.5 logg 9.5 in 
steps of 0.25, -2.0 ^ [Fe/H] ^ 0.0 in steps of 0.25, and 
kms~^ ^ vsini ^ 1000 kms~^ in steps of 50 kms~^. Models 
at intermediate parameter values were constructed by linear 
interpolation on this grid. Pixels with non-zero data quality 
flags set were excluded from our fits, as well as the regions 
around the spectral lines, OI at 1304 A CH at 1335 A and 
the quasi-molecular H feature at ~ 1400 A. 

Before our models could be compared to the observed 
spectrum, they were convolved with a Gaussian filter to the 
spectral resolution appropriate to the COS (< 1900 A) and 
STIS (> 1900 A) data, and linearly interpolated onto the 
observational wavelength scale. In general, our model fits 
did represent the main features of the data quite well, but 
the formal was somewhat high. In order to obtain more 
realistic formal parameter errors, we therefore added an in- 
trinsic dispersion term to the fiux errors in such a way that 
= 1. This dispersion corresponds to about 1 % of the flux 
at 1900 A. In order to explore the systematic uncertainties 
affecting our fits, we tried fitting many different wavelength 
regions, such as the far-UV (COS), near-UV (STIS), regions 
with a high line concentration and also the whole wavelength 
range spanning both the far- and near-UV. 



3.2 Temperature and Metallicity 

A fit to the far UV-region 1268 A - 1900 A yields Tos = 
14200 ± 50 K, logp = 8.2 ± 0.04, [Fe/H] = -1.2 ± 0.05 and 
vsini — 180 ± 20 kms^^. These fit parameters provide a 
good approximation also for the near-UV spectrum out to 
2800 A where the spectrum may begin to show disc features. 
We adopt these parameters as our best-fit parameters. The 
quoted errors are formal errors from the fitting. Table [l] 
presents the best parameters along with their uncertainties. 
These uncertainties are larger than the formal errors, be- 
cause they include our best estimates of the systematic un- 
certainties associated with different choices of fitting win- 
dows. However, there are several reasons for adopting the 
parameters inferred from the fit to the COS far-UV data 
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Figure 3. The top panel shows the mean out-of-eclipse spectrum for the far UV region, where models at three different [Fe/H], 
(green), -1.2 (red) and -2 (blue), are plotted on top. The bottom left panel shows a zoom of the region around the Si II, where the 
coloured lines correspond to the same values of the metallicity as in the top panel. In the bottom right panel, the colours correspond to 
vsini at kms~^ (green), 180 kms~^ (red), 500 kms"^ (blue). 



set as our best-bet estimates. This region includes almost 
all line features, and is also most sensitive to changes in Teff 
and log g. Also, the far UV-region is more likely to represent 
pure light from the WD, while towards redder wavelengths, 
the spectrum could be affected by the disc and the bright 
spot. Finally, fitting only the COS data removes any possi- 
bility that a mis-match in the flux calibration of COS and 
STIS could affect our results. 

The two top panels in Figure [2] show the best-fit model 
in red, along with the mean out-of-eclipse spectrum for 
J1507. The corresponding residuals between model and data 
are plotted underneath each of the COS and STIS spectral 
regions, together with the 1-a error range marked in green. 
The dips in the residuals at (1460, 1530, 1600 and 1670) A, 
are caused by bad pixels in the data and these regions were 
excluded from our fits. As a com parison, a model showing 
the parameters by iLittlefair et a l. (2007), T^s = 11000 K 
and \ogg = 8.5, is plotted in blue. Clearly, this temperature 
is far too low to match the data. 

The top panel in Figure O shows the mean spectrum 
plotted together with models at three different values of 
[Fe/H], (green), -1.2 (red) and -2 (blue). It is clear that 
a model of solar metallicity (green line) does not match the 
data and overestimates both the number and strength of vis- 
ible absorption lines in the spectrum. The left panel shows a 
zoom of the region around the Si II line (the coloured lines 
correspond to the same values of the metallicity as given for 
the top panel) . To the right is a plot showing the dependence 
of the model on t;sinj, where the coloured lines in green, 
red and blue represent rotational velocities of kms~^, 180 
kms~^ and 500 kms~^, respectively. In all panels, the red 
lines show the best fitted parameters. The metallicity and 
rotational velocity affect both the line depths and widths. 
Therefore, as a check, we fixed T^s and logg to our best 
values, and made model fits to smaller wavelength regions 
specifically chosen because they contain strong and/or many 
spectral lines, with v sin i and [Fe / H] as free parameters. We 
find consistent values of vsini and [Fe/H], independent of 
the wavelength region we perform the fits on. 



We compared the metallicity found for J1507 with the 
metallicity distr ibution for Gala c tic ha lo stars at the same 
space velocity. [Patterson et al.l l|2008l ) calculated a space 
velocity for J1507 of 167 kms~^, which can be broken 
into the Galactic velocity compon ents U, V and W, w here 
VIPTW^ =139 kms"^ for J1507. ICarnev et al.ll|l996l ) pre- 
sented a study of the metallicity distribution for the Galac- 
tic halo stars. They find that for stars with velocities in the 
range 100 < VU'^ + < 140, the metallicity distribution 
peaks around -0.7, but has a long tail towards metallicities 
as low as -1.5 (see Figure 6 in their paper). This is in good 
agreement with our metallicity of -1.2. Therefore, we con- 
clude that J1507 is most likely a Galactic halo CV. 



3.3 Possible Extinction Effects 

Up to now, we have assumed that the observed far-UV spec- 
trum of SDSS J1507 is unaffected by interstellar extinction. 
Based on the absence of the well-known 2175 A absorption 
feature in the data, we can set E(B-V) ^ 0.05 as a fairly 
conservative limit on the amount of reddening that may be 
present. In order to test the effect of extinction at this level 
on our conclusions, we carried out additional model fits af- 
ter dereddening the data by E(B-V) = 0.05. We find that 
none of our inferred parameters would change significantly 
if extinction at this level were present. 



3.4 Distance Estimates 

The theoretical relationship between observed flux (F) and 
the Eddington fiux (H) provided by the Synspec models 
can be used to estimate the distance (d) towards J1507. More 
specifically, F = 4nR^^H/ d^, where f?wd is the WD radius. 
This means that the normalisation factor needed to opti- 
mally match a model spectrum to the data is a direct mea- 
sure of R^^/dP . WDs also follow the well-known mass-radius 
relation (with only a weak temperature sensitivity), so the 
surface gravity of the model, g — G M^^/ R^^, uniquely fixes 
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Table 1. Best-fit parameters for J 1507 obtained from 
model fitting to the far-UV region. The errors are de- 
fined as the whole range for which we find good fits, 
irrespective of wavelength region. Formal 1-ct errors 
are given in parenthesis. 





14200 ± 500 (50) K 


logg: 


8.2 ± 0.3 (0.04) 


[Fe/H]: 


-1.2 ± 0.2 (0.05) 


V sin i: 


180 ± 20 (20) kms-l 



J?wd. Thus for a WD model with given Tctt and log the nor- 
malisation factor required to fit the data is a unique measure 
of distance. 

In practice, we estimated the distance by first fitting a 
linear function to the relationship between Mwd and log g in 
Pierre Bergeron's WD cooling model^, at Tea = 14500 K. 
We then use this function to estimate Mwd, and hence J?wd, 
for given log g. Rough values of the WD mass and radius are 
found at Mwd = 0.75 ± 0.15 Mr and i^wd = 0.011 ± 0.002 
Rp). This is co nsistent with both 'Littlefa ir et~all (j2007h and 
[Patterson et a l. (2008). The normalisation constant of the 
model combined with 7?wd , yield the corresponding distance 
estimate. We find a distance towards J1507 of d = 250 ± 50 
pc (the effect of reddening is allowed for in the quoted error). 



4 DISCUSSION AND SUMMARY 

We have obtained HST observations in the UV spectral- 
range of the cataclysmic variable J1507, with the aim of 
measuring its metallicity to determine whether or not the 
system is a member of the Galactic halo. 

By comparing the observed spectrum to synthetic spec- 
tra described by the four parameters, Tes, log g, [Fe/H] and 
vsini, a best fit is found at T^s = 14200 ±500 K. This value 
of TeS is based on the assumption that the WD is the only 
component contributing to the UV flux in J1507. If a hot, 
optically thick boundary layer would be present, it would 
bias our estimate of Tch. However, the single-temperature 
WD models presented here seem to provide a good fit to 
the data, and we expect that the boundary layer in a low 
M-system (such as J1507) would be optically thin and thus 
not contribute significantly to the UV fiux. 

Our best model fit give a T^s that is much higher 
than previous estimates, which implies that the accretion 
rate is higher than previously suggested. More specifically, 
since M oc T*, the 30% increase in the estimated Teff 
corresponds t o an increase in M by almost a factor of 3 
ijTownslev fc B ildstcn 2003). This may have significant im- 
plications for the evolution of this system, including the 
question of whether gravitational ra diation alone is suffi- 
cient to drive this accret ion rate (see iTownslev fc Gansick^ 
I2OO9I : iKnigge et aLll201ll ) . 

At the higher Teff we infer, Fe II and III are no longer 
dominant contributors to the atmospheric UV opacity, mak- 
ing it more difficult than expected to measure the metal- 



^ Coohng models by Pierre Bergeron are found at: 
|http://www.astro.umontreal.ca/ ~bergeron/CoolingModels/ 1 



licity of the system. Nevertheless, model fits to the data 
clearly favour a significantly sub-solar metallicity, [Fe/H] — 
— 1.2 ± 0.2, comparable to the typical metallicity found for 
halo stars at the same hig h space velocity as J1507. 

[Patterson et al.l (|2008l ) found pulsations in J1507, which 
they identified as non-radial pulsations originating from the 
primary WD. Non-accreting CVs show pulsations within 
a well-defined temperature ran ge, the so-called , insta bility 
strip at 10900 K - 12200 K (jGianninas et al.ll2006l ). We 
find that the effective temperature in J1507 is well above 
the instability strip for non-accreting pulsating WDs. How- 
ever, it is already becoming clear that non-radially pulsating 
WDs in CVs are found across a wi der temperature ran ge, of- 
ten towards higher temperatures (|Szkodv et al.|[2010l ). The 
amplitudes of these pulsations are expected to vary with 
wavelength, and a search for the UV counterpart of the op- 
tical WD pulsations will be presented in Uthas et al. (in 
preparation). The fact that non-radial pulsations are present 
in J 1507 is certainly interesting and potentially important 
since a higher metallicity in the outer envelope of the ac- 
cretors might be able to push the instability strip towards 
higher temperatures (| Arras et al.l i2006.). In this context, it 
is interesting to note that J1507 has a very low metallic- 
ity, but is nevertheless managing to pulsate at an effective 
temperature of above 14000 K. 
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